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1. Introduction 



In a recent publication jj] we investigated the role of the negative-energy projections (NEP) of the 
relativistic bound nucleon wave function in the longitudinal (L), transverse (T, TT) and longitudinal- 
transverse (TL) response functions of (e,e'p) processes. Our work was motivated by the success in 
the interpretation of recent (e,e'p) data Q of the relativistic distorted wave impulse approximation 
(RDWIA) M ||, |j| that uses relativistic bound and scattering wave functions (based on the S-V Walecka 
model and Dirac phenomenology). Data on momentum distributions obtained from exclusive 
quasiclastic experiments over extended ranges of missing momenta at low missing energies have been 
successfully analyzed in RDWIA producing reasonable values of the spectroscopic factors (s a ~ 0.7 for 
the S1/2 shell in 208 Pb) consistent at low and high momentum || ||. 

The study carried out in ref. jl| was made in the relativistic plane wave impulse approximation 
(RPWIA) that differs from RDWIA in that the scattered nucleon is described as a plane wave. Both 
RPWIA and RDWIA calculations can also be found in refs. || |l0[ [llj]. The present work focusses on the 
TL response in RPWIA for the reasons discussed below. 

The investigation in ref. |Q] was focused on the relativistic plane-wave impulse approximation (RP- 
WIA) with the purpose of isolating effects due to the relativistic bound nucleon wave function and current. 
We discussed the results for three choices of the nucleon current gauge: Landau or NCCi, Coulomb or 
CCi(O) and Weyl or CCi(3). We also considered two current operators that are called i=l,2 following 
de Forest's notation |l2|. It was found that the role of the NEP of the relativistic bound nucleon wave 
function in RPWIA is to cause a reduction of the differential cross section at low p (p < 300 MeV) and an 
increase at high p (p > 300 MeV) . In RPWIA |l|] , the reduction at low p is generally small and depends 
little on the current operator and gauge choice with one proviso: the Weyl gauge choice may produce 
up to 20% deviations from Landau and Coulomb gauges at low p. However this deviation is also present 
when the NEP of the relativistic bound nucleon wave functions is neglected. It is known to be due to the 
difference between <tcc(o) an d ^cciz) defined in ref. p"2[ , therefore appearing even when non-relativistic 
(n.r.) bound nucleon wave functions are used [jl3| . At larger p values (p > 300 MeV) the size of the NEP 
contributions depends largely on the choice of the current operator and gauge. The NEP contributions 
are minimized for the Landau and Coulomb gauge and for i=2 current operator, while the Weyl gauge 
and the i=l current operator tend to maximize the effect of the NEP contributions, producing large 
differences in some of the response functions. This is particularly the case for the TL response function 
(R TL ) where we found variations of up to a factor of three. 

At present a certain degree of controversy surrounds the TL responses obtained from exclusive 
quasiclastic electron scattering experiments for the least bound orbitals in several nuclei ( 12 C, 16 0, 
208 Pb ): in some cases |l4| large deviations from standard DWIA calculations appear while in oth- 
ers (fj| the data are close to the DWIA calculations. 

Since our study in ref. Q was focused on the P1/2 shell in 16 O, it is interesting to know to what 
extent the above mentioned discrepancy by a factor three remains when knockout from other orbitals is 
considered. Hence in the present work we extend our study to the I753/2 orbital (in 12 C and 16 0), as 
well as to the 1^3/2, Ids/2 orbitals in 40 Ca. The aim here is to see whether the large effect on the R TL 
response found in ref. (l) is a general feature present in other shells, or is particular to the lpi/2 shell 
considered in our recent studies. 

With this goal in mind we restrict ourselves to the RPWIA, as in our previous work. Of course, before 
detailed comparisons can be made with experimental data it is necessary to incorporate the effects of 
final-state interactions (FSI) in a full RDWIA approach. Although FSI are needed for a quantitative 
comparison with experimental data || [|, |j, we expect that, at least for the largest contributions to 
the cross section (i.e., excluding from this statement some of the smaller responses, especially some of 
those that involve polarization degrees of freedom) , the effects of FSI to be basically the same for one- 
nucleon knockout from different bound orbitals, provided the outgoing kinetic energy is the same. At 
the maximum of the TL response that forms the focus of the present work one finds that for instance 
for the pi/2 and p 3 / 2 shells in 16 the FSI effects differ at most by a 4% We therefore expect 

that the conclusions of our study here can be extrapolated to the more realistic situation when RDWIA 



2 



calculations are performed. On the other hand the present study is necessary to identify the dependence 
on the NEP of the bound nucleon wave function for different orbitals. 

The organization of this paper is as follows: in Section 2 we introduce, within the context of the 
Walecka model, the scalar, baryonic and vector densities in coordinate and momentum space. Results for 
various spin-orbit partners are shown and discussed. In Section 3 we study the transverse-longitudinal 
response function for one-nucleon knockout from a selection of relativistic bound orbitals. We compare 
and discuss the results obtained for different current operators and gauges and the role of the negative 
energy projections of the relativistic bound nucleon wave function is studied. Finally in Section 4 we 
present a summary and our conclusions. 



2. Bound nucleon wave functions and densities of spin— orbit part- 
ners 



As is well known, in the Walecka model p, |7j the Dirac equations for bound nucleons in a finite nucleus are 
derived from an interacting relativistic field theory of mesons and baryons by approximating the meson 
field operators by classical fields. For spherically symmetric spin and isospin saturated nuclei treated 
within the context of the Hartree approximation only scalar, S = S(r), and vector, — (V(r),0), 
potentials are to be considered and we write the time-independent Dirac equation with S — V potentials 
as 



jqE — 7 • p — M xp = , 



with 



(1) 



M = M-S 
E = E~V 



(2) 
(3) 



where M is the nucleon mass and E is the total energy eigenvalue (that is, e = M — E > for bound 
states). Here and in what follows the conventions of Bjorken and Drell jl7| are followed. From eq. (Q) 
two Schrodinger-like equations can be derived for the upper and lower bispinor components of ip: 



V t down 



(4) 



with 



1 dA; 
f A + dr 
1 dA. 



a t _ d_ 

r dr 



A+A- 



A- dr V r dr ' 



(5) 
(6) 



A± 

A+A- 



E±M , 

E 2 - M 2 = K 2 



(7) 
(8) 



or 



with 



and 



K 2 = K 2 + 2MU central 



K 2 = E 2 - M 2 



U, 



central 



V 2 - S 2 - 2EV + 2MS 
2M 



(9) 
(10) 

(11) 
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Here we have defined K 2 — E 2 — M 2 < for bound states, and a central potential U ce ntrai that 
depends on the energy E. Eqs.( [sj) and (|j) are Schrodinger-like equations for the upper and lower 
components, respectively, containing in addition to the kinetic energy and central potential (V central) 

terms, a spin-orbit ( — — — ) term and a Darwin (— — ^ ■ — ) term. The last two being generally 

A± or r A± or or 

different for the upper and lower components. 

We recall that in the relativistic case the bound nucleon wave function is labeled by the quantum 
number n that characterizes the eigenvalues of total angular momentum and parity of the four-spinor. 
The relativistic parity operator II = 7 P, with P<fi(r) = <f>(—r), has the eigenvalue II = (— l) K jh and the 
j-eigenvalue is j — \k\ — \- In standard notation tp K is written as 



CM 



ff»(r)C(^) 
ih(r)4>^{r) 



with 



Act 



Act 



The relationship between a, j, £ and £ is: 
For k > 0: 



For k < 0: 



K=j + 



3 + 



£ + 1 , n = {-if 



-(£+!) = -£ , n=(-i) 



(12) 

(13) 
(14) 

(15) 
(16) 



Because in studies of nuclear structure we are more familiar with the non relativistic labeling of 
single-particle states, in practice it is convenient to resort to the £j characterization of the wave functions 
as in the non— relativistic limit, where £ is a good quantum number and the spin-orbit partners are the 
j = £± \ states. To fix ideas and to avoid confusion, in what follows we refer to the "stretched" (j = £+\) 
and "jack-knifed" (j = £ — i) partners as the states with n = —{I + 1) < 0, and K = £ > 0, respectively, 
with £ characterizing the orbital angular momentum of the upper component as in eqs. (p~2| |l6[) . Hence 
the lower components of the spin-orbit partners have £ = £ + 1 and £ = £ — 1 for the stretched and 
jack-knifed £j partners, respectively. 

The general formulas above may be cast in various forms — here we write them as equations for 
upper and lower components with effective central and spin-orbit potentials to bring out the essential 
differences to be found for the radial wave functions that result. Alternative schemes may be used (for 
instance, through the introduction of an effective mass; see, for example, 0), although the present one 
is sufficient for our purposes. The Schrodinger-like equations for the upper components of the spin-orbit 
partners are: 

For the stretched case {j = 1 + 1/2), 



r dr 2 



£{£ + !) 



2MU, 



central 



1 dA + 
A+ dr 



d_ 

dr 



For the jack- knifed case (j — £ - 
1 d 2 1(1+1) 



1/2), 

■f 2Mb central 



1 dA A 
A~~o~r~ 



d_ 

dr 



£+1 



9 k = -K 2 g K (r) 



-K 2 g K {r) 



(17) 



(18) 



Here, in close analogy with the usual non-relativistic limit, the equations for the spin-orbit partners differ 
only in their spin-orbit coupling terms, and these are relatively small. 
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On the other hand the Schrodinger-like equation for the lower component 



r dr 2 



£{£ + !) 



2MU t 



cevitral 



l dA_ 
1Z~dr~ 



j{j+i)-id+i)- 



f K = —Kf K 



(19) 



when applied to the spin-orbit partners reads: 
For the stretched case (j = I + 1/2), 



r dr 2 



(£+!){£ + 2) 



2MU, 



central 



l dA_ 

A- dr 



d_ 

dr 



1 + 2 



for the jack-knifed case (j = £ — 1/2) 



r dr 2 



£{£-!) 



2MU ( 



central 



1 dA_ 

A- dr 



d_ 

dr 



f K = -K 2 f K 



f K = K f K 



(20) 



(21) 



Hence, in contrast to the case of the upper components, the equations for the lower components 
differ not only in the spin-orbit coupling term but also in the centrifugal barrier. This produces stronger 
differences on the radial dependence of the lower components than on the radial dependence of the upper 
components. This is to say, if one compares for instance the radial functions f K for the pi/2 and p$/2 
cases one expects to see larger differences than if one compares the g K for those same shells. 

Later we shall see the consequences of these differences in upper and lower components when examining 
the RPWIA results in detail; here it should be noted that in general the lower components are small and 
the radial dependence of both scalar and baryonic densities is mostly dominated by the upper components: 



Ps(r) 



Pb(t) 



4tt 



2j 



4-7T 



2.7 



-J2^ + ^n= [gUr)+f K (r)} , 



with 



pB{r)r 2 dr = 1 



(22) 
(23) 

(24) 



In the cases considered here the contribution of f K to the normalization is less than 3%. Nevertheless 
it is important to realize that for spin-orbit partners the lower components can differ significantly, while 
the differences in the upper components tend to be mild. This is illustrated in fig. 1 where we plot 
the sum and difference of baryonic and scalar densities (in fm -3 ) for the spin-orbit partners P1/2, P3/2 
in 16 and d 3 / 2 , db/2 m 40 Ca. The 3si/2 orbital in 208 Pb is also plotted for comparison. The present 
calculations follow similar lines to those of ref . |lj ; in particular we use the parameters of ref. |l8| and the 
TIMORA code |l^] . Clearly the differences of f K due to the different centrifugal barriers are maximal in 
the Pi/2~P3/2 spin-orbit partners because the lower component for the pi/2 partner is an s-wave while for 
the P3/2 is a d-wave. The d^/2~d^/2 spin-orbit partners have correspondingly a p-wave and an f-wave as 
lower components. It is also interesting to see that for the 3si/2 case the lower component looks negligibly 
small and may be expected to play a minor role. 

We note that the elementary three-vector current density 



v = tp ~/xp K = v(r)i) 
is proportional to the product of g K and / K 

v(r) = 2g K (r)f K (r) . 



(25) 



(26) 



This quantity is also plotted in fig. 2 for the same orbitals shown in fig. 1. 

The larger differences seen for the lower components reflect a genuine relativistic effect that cannot be 
accounted for when non-relativistic bound nucleon wave functions are used. This naturally motivates an 
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investigation as to whether or not this relativistic effect has observable consequences. Electron scattering 
at intermediate energies, involving large momentum transfers, is the ideal place to look for such relativistic 
effects. 

In elastic electron scattering one is sensitive to the Fourier transforms of Pb{t) and of v(r) and 
relativistic calculations have been presented for a few elastic magnetic form factors [p0| . A thorough 
discussion of relativistic calculations of transverse form factors for elastic electron scattering could provide 
the motivation for future work — however, following our recent studies Q the focus of this paper is 
instead placed on the longitudinal-transverse response functions (R TL ) for exclusive quasielastic electron 
scattering (e,e'p). This particular response function depends on linear combinations of scalar, baryonic 
and vector densities in momentum space 

ps(p) = j^j;?:(p«(p)=d)-^) 1 (27) 

Arrr 

Pb(p) = ^— l Y.^n{pW m {p)=9l{p) + !l{ P ), (28) 

v(p) = -^-^^( P )p.7C(p) = 2nfc(p)/^), (29) 

with ip^lp) the Fourier transform of ip^r). The radial functions in momentum space are given by (see, 
for instance, |EJ) 



9x(p) = \ ~ / r drg K (r)ji(pr) , (30) 



UP) = y-J Q r J drf K (r)jj(pr) . (31) 

Specifically, the relevant quantities in p-space that enter directly in the calculations of response 
functions for (e,e'p) processes are the amplitudes a K and (3 K of the positive and negative energy projections 
of the bound nucleon wave function in the expansion: 

r m {p) = E [& K (p)n(p,s)(s\i^(p)) +0 K (p)v( P ,s)(s\i i ct>- K (p))] , (32) 

S 

with 



a-w - V^f(if¥"- (p) -ra / - (f,) ) ■ (34) 

where u(p, s) and v(p, s) are the positive and negative energy solutions of the free Dirac equation with 
E = \J M 2 +p 2 . The absolute values of the positive and negative energy projection amplitudes (in fm 3 ) 
are plotted in fig. 3 for the spin-orbit partners Pi/2 - i>3/2 and d 3 / 2 -d 5 /2- Note the different scales of the 
positive and negative energy projections. Note also that in the limit of free nucleons f3 K {p) = 0, and 
therefore \(3 K \ 2 measures the dynamical enhancement of the lower components. 

A much larger dynamical enhancement is observed in fig. 3 for the jack-knifed states than for the 
stretched states. This was expected from our previous discusion on the Schrodinger-like equations for 
the upper and lower components. As a matter of fact, much can be learnt by looking at the Schrodinger- 
like equations that follow from the Dirac eq. (|l|). For instance, as first pointed ouy by Ginocchio 
the Schrodinger-like equation for the lower component is the same for relativistic states having the 
same ?-value but £- values differing by two units (as for instance the s 1 / 2 -d 3 ^2 or P 3/2~ /s/2>-" pairs) 



when the spin-orbit coupling term in eq. (19) is negligible. This has been identified as the origin of 
pseudospin symmetry in nuclear spectra by Ginocchio j2^] . Another interesting limiting case can be that 
of quasidegeneracy of levels with equal j and oposite parity that is also present in nuclear spectra [B2j . 
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3. Longitudinal— transverse response function 



In this section we compare results on the R TL response functions for one-nucleon knockout from various 
relativistic bound orbitals. The details of the calculations are as in ref. |ij and we do not repeat them here; 
we use what was called kinematics I in that work, namely quasielastic conditions with q — 500 MeV/c 
and u> = 131.56 MeV. We recall that, using projection techniques, the R TL response can be separated 
into a contribution from the positive energy projection (i?p L ), a contribution from the negative energy 
projection (Rjf), and a crossed term {R^): 

RF^ = i?p^ + R^f ~l~ > (35) 

with 

R T P L = 1l T u tN uu {p) , (36) 

R T N L = TZ T ^N vv (p) , (37) 

R T C L = 1Z^N uv (p) . (38) 

In eqs. (p6|-|38|) the dependence on the nuclear structure is factorized. This dependence is contained in 
the bound nucleon momentum distributions 

N uu (p) = {& k { P )) 2 /Att , (39) 

N VV ( P ) = (4(p)) 2 /4tt , (40) 

N uv (p) = -2& K {p)P K {p)/4ir . (41) 

In contrast, TZ^, TZ^ and H^v depend only on the current operator and gauge and are independent of 
the nuclear structure. 

When the nuclear structure is treated in the non-relativistic limit the nucleon current operator is 
expanded in the basis of positive energy free Dirac spinors (see for instance ref. p3fl ) and the single- 
nucleon response function is 

k tl = n T J: , (42) 

corresponding to electron scattering from a free (relativistic) nucleon, and 

R TLn^ R TL (43) 

With relativistic bound nucleon wave functions one has additional single-nucleon responses TZ^v and TZ^ 
that do not appear in the scattering of electrons from free nucleons (or antinucleons). The appearance 
of these terms is due to the nonzero overlap of the bound nucleon wave function with the Dirac sea that 
results in finite values of the N uv (p) and N vv (p) functions. 

In ref. [pj we showed that the single-nucleon responses and TZ^ change much more with the 
changes of the current operator and gauges than does the TZ 1 ^ response. Whereas the single-nucleon 
response TZ^ has often been studied Q , the TZ^ and 7£jjf were for the first time identified and studied 
in ref. 0. The fact that these responses are more strongly dependent on the choice of the current 
operator and gauge, results in a larger theoretical ambiguity in the calculated total response R TL when 
using relativistic bound nucleon wave functions. 

As we showed in the previous section, the overlap with the Dirac sea of the bound nucleon wave 
function depends strongly on the particular tj orbital under study. Therefore a different sensitivity of 
the R TL response to the contributions from the negative energy projections can be expected for different 
orbitals. In particular we show in figs. 4 and 5 that this sensitivity is quite different for the spin-orbit 
partners pi/ 2 and p 3 / 2 or d 3 / 2 and d 5 / 2 . 

The longitudinal transverse response functions for one-nucleon knockout from the p\/ 2 and p 3 / 2 
orbitals in 16 are shown in the top and bottom panels of fig. 4. For each shell results corresponding 
to the CC1 current operator and to the three different gauges (Landau, Coulomb and Weyl) are shown 
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on the left panels, while the right hand panels contain the corresponding results for the CC2 current 
operator. 

In these figures one can see that for fixed current operator and shell the Landau and Coulomb gauges 
give practically the same results, while the Weyl gauge tends to produce important deviations. The 
importance of these deviations is much higher for the than for the P3/2 shell. This is also the case 
when changing the current operator. The following ratios are found between the maximum values of the 
R TL response of the p 1 / 2 shell: if we fix the gauge to be the Coulomb gauge, the ratio between the results 
for the CC1 and CC2 current operators is 

R TL (CC2(0)) ~ ' ( ' 

and the ratio is larger in the Weyl gauge 

R TL (CC1(3)) 

R TL (CC2('S)) ' 1 ' 

Similar ratios are found when we fix the current operator and compare the results of different gauges 

R TL (CC2(3)) ^ R TL (CC1(3)) ^ 

R TL {CC2(0)) ~ ' R TL {CC1(0)) ' [ ' 

leading to the large difference, by as much as a factor of 3, between the CC1(3) and the CC2(0) results, 
already mentioned in ref . jlj . 

On the contrary, for the spin-orbit partner P3/2 the results between different choices are not that 
large so the above ratios reach at most a value of 1.2 (that between CC1(3) and CC2(0) results) and in 
most cases are less than a 10%. In particular for the P3/2 case the dependence on the gauge is almost not 
visible with the CC2 current operator. 

We note that R TL contains the product of transverse and longitudinal currents. The transverse 
current changes only with the current operator (i.e., CC1 or CC2) but it is independent on the gauge. 
Therefore, for the cases considered, the strong gauge dependence for the p 1 / 2 case must come from the 
longitudinal current. A similar comparison is made in fig. 5 for the <i 3 / 2 and d 5 / 2 spin-orbit partners. 
It is also the case that the dependence on the current operator and gauge is more pronounced for the 
jack-knifed than for the stretched cases. 

As seen in figs. 6 and 7 when we consider only the contribution from the positive energy projection 
(i.e., the truncation in eq. (^j)) the sensitivity of the R TL response to the choice of the current operator 
and gauge is similar for the spin-orbit partners P\/2~Pz/2 and d3/2 _ ^5/2- The larger sensitivity obtained 
for the P1/2 and d 3 / 2 partners is due to their larger negative energy projections at low p (p < 200 MeV) 
that carries larger deviations particularly in the R^ L term. This is in turn an effect of the Schrodinger- 
like equations for the lower components of the spin-orbit partners. 



4. Final remarks 

The main conclusion of this work is that the large deviation Jl|] in R TL predictions produced by different 
current operators and gauges takes place for the jack-knifed states (pi/ 2 , d 3 / 2 ,...), but not for the stretched 
states (p3/2j ^5/2v)- This genuine relativistic effect stems from the dynamical enhancement of the lower 
components, and can be traced back to the Schrodinger-like equations for the upper and lower components 
of the relativistic bound nucleon wave function. 

We have shown that this effect can be understood from the quite different behaviour of the dynamical 
enhancement function /3 K of the stretched and jack-knifed states. Indeed, from eq. (|34|) and fig. 3 one can 
see that from the jack-knifed states (i.e. Pi/ 2 , d 3 / 2 , ■ ■ ■) the amplitudes of the negative energy projections 
are much larger than those for the stretched states (p3/ 2 , ^5/2, ■ ■ ■)■ O ne has to realize that in eq. ( |34|) 



Pk(p) goes like the difference of two functions of p with the same number of nodes (n + £ + 1) for the 
stretched states, while for the jack-knifed ones, it goes like the difference of a function with n + £ + 1 
nodes and a function with n + £—l nodes (n = 1 in the cases considered here) . This is due to the different 
quantum number £ of the lower components in stretched and jack-knifed states and can be traced back 
to the Schrodinger-like equations (|2C|)-(|2l|). This is in contrast to the amplitudes of the positive energy 
projections (oe K ) that in all cases are strongly dominated by the upper component (g K ) and are therefore 
similar for stretched and jack- knifed states. Note that g K (r) for stretched and jack-knifed states (with 
equal n and £) differ only in their r.m.s. radii. 

Since the non-relativistic limit corresponds to f3 K (p) = 0, as a corolary we may also conclude that 
comparing to the non-relativistic limit, the responses for the stretched states will be closer than those 
for the jack-knifed states to their respective non-relativistic limits. 

The choice of current operator and gauge is an important issue in quasielastic electron scattering 
since theoretical results on differential cross-sections and response functions depend substantially on 
those choices (see for instance refs. jy|, ^3|, Q and refs. therein). This tends to be particularly so when 
relativistic bound and/or scattering wave functions are used because the contributions from the negative 
energy projections depend more strongly on the current operator choice than do the contributions from 
the positive energy projections. The relativistic approach provides a tool to study further reaching 
consequences of the various choices of the current operators and gauges that can not be explored if one 
attaches to the non-relativistic description of the nuclear structure and wave functions where one can 
not go beyond the truncated expression in eq. (|43|). 

In the case of stretched states the differences in the R TL responses with different current operators and 
gauges are of the same order (and have the same origin) as those found with non-relativistic bound nucleon 
wave functions. The theoretical uncertainty for these states is at most a 20% (this is the largest ratio 
between CC1(3) and CC2(0) calculations), which is small compared to the large theoretical uncertainty 
of up to 300% found for the jack-knifed states. Although our calculations here do not contain final-state 
interactions, a similar behaviour is to be expected when FSI are included. Since FSI and short-range 
initial-state correlations are expected to play a similar role when calculating the R TL response of different 
spin-orbit partners, the effect found here may provide a way to test different choices of current operators 
and gauges. 

Since the results for jack-knifed states amplify differences between models, while those for the stretched 
states are more "model independent" , the combined analysis of data on spin-orbit partners can be used to 
elucidate between models. On the other hand, if one wants to determine with minimal model dependence 
a certain observable, one should focus on measurements on stretched states. 
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Figure captions 



Figure 1: Semi-sum (left panels) and semi-difference (right panels) of the baryonic and scalar densities. 

Results are shown for the spin-orbit partners: lpi/2 (thick-dotted line), lp 3 /2 (thin-dotted line), 
and ld 3 / 2 (thick-dashed), 1^5/2 (thin-dashed). For comparison we also present results for the shell 
3si/2 (solid line). 



Figure 2: Same as fig. 1 for the vector current density (eq. (26)). 



Figure 3: Amplitudes |a K (p)| (top panel) and |/3 K (p)| (bottom panel) for the shells: lpi/2, lP3/2> 1^3/2 an( i 
1g?5/ 2 - The labeling of the various curves is the same as in fig. 1. 

Figure 4: Interference longitudinal-transverse response function R TL (in units of fm 3 ) for the shells: lpi/% 
(top panels) and IP3/2 (bottom panels) in 16 O. Panels on the left hand side correspond to results 
for the CC1 current operator and the gauges: Landau (solid line), Coulomb (dotted) and Weyl 
(dashed). Right hand panels show the results for the CC2 current operator and the same labeling. 

Figure 5: Same as fig. 4 for the spin-orbit partners: 1g?3/2 and 1^5/2- 

Figure 6: Components of the transverse-longitudinal response, i?p L , R^ L and Rjf (see eqs. (|3^ |3§|)). Results 
are shown for the spin-orbit partners: lpi/2 (top panel) and (bottom panel). Thick lines 

correspond to gauges that use the CC2 current operator and thin lines correspond to the CC1 
operator. Gauges considered are: Landau (solid lines), Coulomb (dotted lines) and Weyl (dashed 
lines). 

Figure 7: Same as fig. 6 for the shells: ld 3 / 2 (top panel) and ld 5 / 2 (lower panel). 
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Figure 1: Semi-sum (left panels) and semi-difference (right panels) of the baryonic and scalar densities. 
Results are shown for the spin-orbit partners: lpi/2 (thick-dotted line), 1^3/2 (thin-dotted line), and 
1^3/2 (thick-dashed), ld 5 / 2 (thin-dashed). For comparison we also present results for the shell 3si/2 
(solid line). 
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Figure 3: Amplitudes |& re (p)| (top panel) and \j3 K (p)\ (bottom panel) for the shells: lpi/2, IP3/2 
and Ids/2- The labeling of the various curves is the same as in fig. 1. 
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Figure 4: Interference longitudinal-transverse response function R TL (in units of fm 3 ) for the shells lp\/2 
(top panels) and (bottom panels) in 16 0. Panels on the left hand side correspond to results for the 

CC1 current operator and the gauges: Landau (solid line), Coulomb (dotted) and Weyl (dashed). Right 
hand panels show the results for the CC2 current operator and the same labeling. 
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Figure 5: Same as fig. 4 for the spin-orbit partners: ld 3 / 2 and ld 5 / 2 - 
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Figure 6: Components of the transverse-longitudinal response, i?p L , Rc L and Rjj L (see eqs. (||^ |3S|)). 
Results are shown for the spin-orbit partners: lp±/2 (top panel) and (bottom panel). Thick lines 

correspond to gauges that use the CC2 current operator and thin lines correspond to the CC1 operator. 
Gauges considered are: Landau (solid lines), Coulomb (dotted lines) and Weyl (dashed lines). 
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